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Renormalization of the spin-wave spectrum is discussed in a cubic ferromagnet with dipolar 
forces at Tc 3> T > 0. First 1/S-corrections are considered in detail to the bare spectrum 
ek = \/ Dk^{Dk^ + Sloq sin^ 6^^), where D is the spin-wave stiffness, is the angle between k 
and the magnetization and ujq is the characteristic dipolar energy. In accordance with previous 
results we obtain the thermal renormalization of constants D and (jJq in the expression for the bare 
spectrum. Besides, a number of previously unknown features are revealed. We observe terms which 
depend on azimuthal angle of the momentum k. They are proportional to the square invariant 
\kx — k^l/k in the case of simple cubic lattice and to the triangular invariant Im([fc2, -I- ikyf')/k'^ in 
face- and body-centered cubic lattices, where we assume that z-axis is directed along magnetiza- 
tion. It is obtained an isotropic term proportional to k which makes the spectrum linear rather than 
quadratic when sin^k = and k <^ loq/D. In particular a spin- wave gap proportional to sin^k is 
observed. Essentially, thermal contribution from the Hartree-Fock diagram to the isotropic correc- 
tion as well as to the spin-wave gap are proportional to the demagnetizing factor in the direction of 
domain magnetization. This nontrivial behavior is attributed to the long-range nature of the dipolar 
interaction. It is shown that the gap screens infrared singularities of the first l/S'-corrections to 
the spin-wave stiffness and longitudinal dynamical spin susceptibility (LDSS) obtained before. We 
demonstrate that higher order 1/S-corrections to these quantities are small at T <^ cjq. However 
the analysis of the entire perturbation series is still required to derive the spectrum and LDSS when 
r ^ liJo- Such an analysis is out of the scope of the present paper. 

PACS numbers: 75.30.Ds, 75.10.Jm, 75.10.Dg 

I. INTRODUCTION 

Despite enormous amount of studies on the theory of ferromagnets with magnetic dipolar interaction there are 
still some unresolved problems. Infrared divergences of the longitudinal dynamical spin susceptibility^'^ and infrared 
divergences of l/S'-corrections to the spin- wave stiffness^ are among them. 

It was obtained by Toperverg and Yashenkin^ that dipolar forces lead to strong long-wavelength fluctuations 
manifesting themselves in infrared divergence of the first perturbation corrections to the uniform longitudinal spin 
susceptibility: X|| {(^ ~^ 0) ^ iT/uj. It is pointed out in Refi that this singularity is nonphysical one since it leads to a 
nonzero value of absorption function cx ujliaxiuj) at a; = and thus the sample would be heated by a dc field. Then 
one could expect great renormalization of X||(a;, k) after taking into account higher order perturbation corrections at 
small energy and momentum. Thus, experiment on nearly isotropic ferromagnet CdGr2Se4 shown weaker power-law 
dependence*? X\\{^ ^ 0) ^ {i/uj}'^-'^^. The analysis of the entire perturbation series for X||(w, k) has not been carried 
out yet. 

First 1 / 5-corrections to the spin- wave spectrum have been discussed by Rahman and Mills in Refi^. They found an 
infrared divergent contribution to the spin-wave stiffness. It stems from three-magnons terms in the Hamiltonian. As 
they did not obtain any term remaining finite as /c — > 0, one could expect a great renormalization of the spectrum at 
small momenta. However Rahman and Mills used an approximation discussed below in detail to simplify intermediate 
_ ' calculations. It remains unclear if there exists a constant term as fc — s- outside this approximation. 
k> These problems of infrared singularities are related to the fact that a spin-wave gap has not been obtained yet in 

5^ magnon spectrum. It is well known that within the linear spin-wave approximation the spectrum has the following 
form at small fct^ 



Ek = y [Dk"^ + g^iHi')) {Dk^ + g^iHd) + Sluq sin^ 0^) , (1) 

where D is the spin-wave stiffness, iJ^'-* is the intrinsic magnetic field, is the angle between k and the magnetization 
and 

.0 = (2) 

is the characteristic dipolar energy, where ■^o is the unit cell volume. It is well established that a finite-size sample of 
a ferromagnet breaks up into domains at zero external magnetic field so as the intrinsic magnetic field is zero.^ It is 



2 



seen from Eq. |^ that the spectrum has no gap in this case. The intrinsic magnetic field is also zero in a unidomain 
sample that is infinite along the direction of magnetization. 

At the same time there are some considerations in favor of existence of the spin-wave gap in a ferromagnet with 
dipolar forces. First, due to its long-range nature and symmetry the dipolar interaction violates the well-known 
Goldstone theorem and can lead to the gap. Second, it was obtained in Ref.^ that pscudodipolar interaction which 
has the same symmetry as the long-range dipolar interaction leads to an energy gap in antiferromagnets. Third, it 
is well-known that within the first order of 1/5 dipolar and pscudodipolar forces lead to anisotropic corrections to 
the total energy of a ferromagnet A^iiS As a result, directions along edges of the cube are energetically favorable in 
a simple cubic lattice (SCL) whereas the magnetization should be parallel to a body diagonal of the cube in the 
face-centered cubic lattice (FCCL) and the body-centered cubic lattice (BCCL)j ^i^i-^°i^^ It has been pointed out by 
Keffer^iii that the anisotropic terms in the total energy of a ferromagnet should be accompanied with an "energy 
shift" in the spin-wave spectrum. However the particular calculations have not been realized yet thoroughly because 
of their cumbersomeness. 

In the present paper we carry out a comprehensive analysis of the real part of first 1 / S'-corrections to the spectrum 
of a ferromagnet with a cubic lattice. We use a technique based on introduction of anomalous Green's functions 
rather than on conventional Bogolyubov's transformation. This technique proved to be more convenient making 
intermediate calculations more compact 1^*^^*^^*^^ In accordance with Ref.'^ we obtain the thermal renormalization of 
constants D and loq in Eq. (^). Besides, we reveal a number of previously unknown features. In particular it is 
shown that there is a gap in the spectrum proportional to sin^^- This gap screens singularities of 1 /S'-corrections 
to the spin- wave stiffness and longitudinal spin susceptibility. We show that higher order l/S-corrections to these 
quantities are small at T <C Wq- However the analysis of the entire perturbation series is still required when T ^ wo 
which is out of the scope of the present paper. In this case one could expect large renormalization of the longitudinal 
dynamical susceptibility as the nontrivial power-law dependence of X|| (lo) discussed above was obtained experimentally 
in CdCr2Se4 at T » wqi^ 

Then, we observe corrections to the spectrum which depend on azimuthal angle of the momentum k. They are 
proportional to the square invariant |fc^ — ky\/k in the case of SCL and to the triangular invariant Im([fc2; -f iky]'^)/k'^ 
for BCCL and FCCL, where we assume that z-axis is directed along magnetization. This difference is accounted for 
the fact that due to the dipolar anisotropy magnetization is parallel to an edge of a cube (tetrad axis) in SCL and to a 
body diagonal of a cube (triad axis) in BCCL and FCCL. It is also shown that there is an isotropic term proportional 
to k that plays the predominant part when sin^?k = and k <C loq/D: renormalized spectrum is linear in k rather 
than quadratic [cf . Eq. J^l] • We demonstrate a relation between the dipolar anisotropy and the gap and isotropic term 
in the spectrum. Both the gap and the isotropic term are proportional to sums over momenta containing anomalous 
Green's functions that were omitted in Ref.'^. 

Remarkably, thermal contribution from the Hartree-Fock diagram to the isotropic correction as well as to the 
spin-wave gap are proportional to the demagnetizing factor in the direction of the domain magnetization. Then a 
spin- wave excitation with wavelength much smaller than the size of the domain appears to "feel" the domain shape. 
This nontrivial behavior is attributed to the long-range nature of the dipolar interaction. 

The rest of this paper is organized as follows. The Hamiltonian transformation and the technique are discussed in 
Sec.m First l/S-corrections to the spin- wave spectrum are considered in Sec. lIIII In Sec.|^we discuss the screening of 
infrared singularities of longitudinal spin fluctuations and further l/5'-corrections to the spin- wave stiffness. Possibility 
of experimental verification of the results obtained is also considered in Sec. lIVI Sec. IVl contains our conclusions. One 
appendix is included with details of calculations. 



II. HAMILTONIAN TRANSFORMATION 



The Hamiltonian of a ferromagnet with dipolar interaction in magnetic field H has the well known form: 




(3) 



Im 



(4) 



We assume that H = in the case of the unidomain sample that is infinite in the direction of magnetization. For 
multidomain sample the Hamiltonian Q describes one domain and H is the magnetic field produced by all other 
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domains. Taking the Fourier transformation we have from Eq. 



(5) 



where Jk — >^im exp(jkRim), Q'^^ = X)/ ^^PC^'^-f^'m) ^^"^ ^ ^^c number of spins in the lattice. Dipolar 
tensor Q'^^ possesses the well known properties which are the same for all types of cubic lattices and independent of 
the orientation of the coordinate frame relative to cubic axes:^'^^ 



Wo 



1 



if y < A: < 1, 

1j 



(6) 
(7) 



where we set the lattice spacing equal to unity, L is the characteristic length of the domain, A/'q are demagnetizing 
factors and loq is the characteristic dipolar energy given by Eq. (j^J . 
After Dyson-Maleev transformation 



( flk + a_k 



25 



"^k 



'-k 



2S 



(8) 



Hamiltonian Q has the form H = + where Eq is the ground state energy and Hi denote terms containing 

products of i operators a and . As we intend to calculate corrections of the first order of 1/5, we need terms up to 
Ha- Hi = because it contains only Q^^ with (3. For the rest necessary terms one has: 



H2 
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t ^Ir t t 

-Ek^flk + — aka_k ^ " " 
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as, 



ki+k2+k3+k4=0 

where we drop index k in Eqs. pU|) and and 
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(9) 
(10) 

(11) 



Bk 



5(Jo - Jk) - 1 ( g 



XX 

k 



- gi.i{4TTgfiSJVz + H) 



+ ^ sin^ 6'k - 9^i{H + 47r5^5A/; 



(Qr - 2^Q;^^ - Qf ) ^ sin^ O^e'^^^^, 



(12) 
(13) 



where the expressions after « are approximate values of corresponding quantities at fc ^ 1. In the multidomain 
sample the expression in the brackets of the last term in Eq. (|12f) is the intrinsic magnetic field which is zero. In 
the unidomain sample that is infinite in the direction of magnetization we have = 0, H = and the last term in 
Eq. p2|) also vanishes. In the spin- wave approximation the magnon spectrum is 



Ek 



El 



(14) 



Using Eqs. I©, ((T^ and ifT^ we obtain from Eq. at L ^ <C fc <C 1 the well known result Q). 

To perform the calculations it is convenient to introduce the following retarded Green's functions: G(ci;,k) 



(ak,aj^)i^, -F(w,k) = (ak,a_k)(^, G'(w,k) = 



-k' 



i-k)c. = G*i~io,~k) and Ft(a;,k) = (aLk:«l)- = F*{-u;,-k). 



We have two sets of Dyson equations for them. One of these sets has the form: 

G(w,k) = G^°\uj,k)+G''"\iJ,k)^Lj,k)Giiu,k)+G<'"^i^,k)[Bl + U{ij,k)]F^Lj,k), 
F\iu,k)=G (w,k)I](aj, k)F^Lj, k) + G^°^ {uj, k)[Bk + nt(cj, k)]G{uj, k). 



(15) 
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where G'-°^ {cj, k) = {uj-Ek + iS)''^ is the bare Green's function and E, S, 11 and IT are the self-energy parts. Solving 
Eqs. p5(l one obtains: 



G(u;,k) 



V{LU,k) ' 



= PRk^ ' ^^^^ 

. i3k + nt(c.,k) 

^ = P(c.,k) ' 

where 

V{u, k) = {uj + i5f -el- Ek{i: + E) + B^n + B^tf + + - S) + nnt - (17) 

and ek is given by Eq. (|14|l . The last two terms in Eq. (|17|l give corrections of at least second order oil/S and will not 
be considered below. It is convenient for the following to introduce a new quantity f2(a;,k) describing the spin-wave 
spectrum renormalization [cf. Eq. H17|l ]: 

V{uj,k) = (uj + iS)^ -el~n{Lu,k). (18) 

We calculate ft{uj,k) in the next section. 

III. RENORMALIZATION OF THE SPIN- WAVE SPECTRUM 

Corrections to the spin- wave spectrum to be obtained are proportional to sums over momenta containing components 
of the dipolar tensor Q"^^ ■ In some of these sums summation over small momenta is important and one can use 
expressions © and Q for Q^^ ■ As properties of the dipolar tensor 10 and {T)) are the same for all types of cubic 
lattices and independent of orientation of the coordinate frame relative to the cubic axes, such sums do not depend 
on the direction of the quantized axis and the lattice type. Meantime there are sums in which summation over large 
momenta is essential and which, consequently, depend on the direction of quantized axis and the lattice type. Thus, 
one should bear in mind what is the direction of magnetization in the ground state. 

Therefore the well known fact should be taken into account that dipolar and pseudodipolar interactions lead the 
energy of a ferromagnet to be dependent on the direction of quantized axis;^'^ It has been shown in Ref.^ that the 
first 1 / ^'-correction to the energy Eq gives rise to such anisotropy. In our notation this correction has the form: 

AE=in,)=j:'-^^-j:^-^, (19) 

k k 

where the isotropic part of the energy is omitted in the right part of this expression. After direct calculations one 
obtains: 

^ = C^hh'y + ihl + (20) 

r, SJ ^{Q^r^Q^f^Hq^ 

^ - ' ('^^ 

where 7i are direction cosines of the magnetization and components of the dipolar tensor in Eq. \'2,1\ are taken relative 
to cubic axes. The constant C can be calculated numerically using the technique of dipolar sums computation (see, 
e.g., Ref.^^ and references therein). We obtain in accordance with Refs.^'^ that it is positive for SCL (C « 0.012) 
and negative for FCCL (C ~ —0.005) and BCCL (C « —0.04). Then, for SCL an edge of the cube is easy direction 
whereas a body diagonal of the cube is easy direction for FCCL and BCCL. 

SCL is considered in detail in the next three subsections (jlll AlUlI Bl and Fill Cp . Calculations and results for FCCL 
and BCCL resemble those for SCL. They are discussed in Sec. IIII Dl 

We study now separately bubble diagrams shown in Fig.^a) a-nd b), and the loop diagram presented in Fig. ^c). 
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A. Bubble diagrams 



Let us begin with the diagram shown in Fig. ^a). It appears from three-magnons terms ()10|) and gives zero. To 
demonstrate this we make all possible couplings of two operators a and in Eq. H1U|) : 



E 

q 



\<^0[Qq +«Qq)+ao(Qq " ) f 



(1 + 2Nq) - ao 



(1 + 2iVq) 



(22) 



ret 



^xz 



where TVq — [e!^^^^ — 1) ^ is the Plank's function. Expression (|22|l is zero because Q 
Eq = Eq: and Bq = Sq/, where q = [qx,qy,qz) and q' = (-^a;, -g^,, g^). 

The Hartree-Fock diagram presented in Fig.^b) comes from 7i4-terms given by Eq. After simple calculations 

we obtain for contribution to f2(Li;,k) from the four-magnons terms: 



f]('*)(w,k) 



Eq{l + 2Nq 



q 

|i3qP(l + 2jVq) 

Se„ 



( Jk — Jo + "^q 



isk 



Eqil 



2Nq) - eq 



q 



Jk-q) + Q^" - Q 



+ Q 



k-q 



Re(B*Bk)(l 



27Vq) 



(23) 



At zero temperature iVq = in Eq. H23|l . In this case the spectrum is renormalized by quantum fluctuations only. 
As the temperature increases, corrections from terms in Eq. (|23l) containing A'q become larger. They exceed terms 
in Eq. (|23|) not containing A^q above a certain temperature. We find below that this temperature is of the order of 
(Duq)^^^. Regimes of T < {Duj^y/^ and T > {Dlj^)^^^ are called below as quantum and thermal ones, respectively. 
It is convenient to discuss them separately. 

T ^ {Duj'^Y^^. To our knowledge corrections to the spin- wave spectrum in this regime have not been considered 
before. Taking A'q = in Eq. H23() one notes that summation over large momenta is essential and {Eq — eq)/eq ~ 



\Bq^' 



As a result, we obtain in the leading order of ujq 

SJ v-^ Q : 



-Dfc^-^sin* 6'k cos^ 20k^ 



q 



Sea 



■ loq sm 



2ea ' 



X = 



^yy 
'q 



q 



4eo 



cos Qx 



0.013, 



(24) 



(25) 



where (/)k is the azimuthal angle of k and X is calculated numerically. The first and the last two terms in Eq. H24|) 
stem, respectively, from the last and the first terms in Eq. (|23|) . Dependence of D,^*^{uj, k) on k differ drastically from 
that of the bare spectrum The first and the second terms in Eq. H24(l give corrections to the spin- wave stiffness. 
It is seen that the first term depends on azimuthal angle (/)k. Such a dependence is allowed by the symmetry because 
the first term can be expressed via the quadratic invariant: sin^ 6'k cos^ 20k — (k^ — ky)^ /k^. Essentially, the second 
term does not depend on direction of k. Thus, it plays the predominant part when sin 0k = and k ujq/ D. In this 
region Q^'^\LL),'k.) is proportional to k^, whereas cx fc^, and the renormalized spectrum appears to be linear rather 
than quadratic. The last term remains finite as k —^ and sin 0k ^ and gives a contribution to the spin- wave gap. 

We stress that in all sums in Eq. (|24|l summation over large q is essential and it is important to take into account 
orientation of quantized axis relative to cubic axes. 

{DujQy/^. Calculation of the terms in Eq. 1)23(1 containing A^q leads us to the following expression in the 
leading order of -Oq- 



17(4) (c^,k) = -{Dky 



-WiT) - Dk'^ujo sin^ ek[WiT) + ViT)] 



W{T) = 



V(T) = 




2uj^AfzV{T) + ^Y. 

q 

5/2 



Sea 



(1 + 2iVq) 



(26) 



(27) 



(28) 
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where w = 0.01 and v = 0.06. The first term in Eq. (|23|) gives the largest thermal contribution to ri('*)(a;, k). Let us 
make sure now that the thermal corrections given by Eq. (|26|l are larger than the quantum ones H24|l at T ^ (DluqY^^. 
Indeed, the first two terms in Eq. H26|l are greater than the first term in Eq. 124|l at such temperatures. Let us compare 
the third term in Eq. 1)26(1 with the second and the third terms in Eq. H24|) . If Afz ^ the first term in the second 
brackets in Eq. are larger than the second one and the last two terms in Eq. HH) when T > {Doj^Y^^. If 
Afz = only the second term remains in the second brackets in Eq. H26|l . It is proportional to loq and TD{ujo/ D)"^^^ 
at T <C V^qD and T ^ ^/ujqD, respectively. We demonstrate below that there is a correction from the loop diagram 
that cancel the T-dependent contribution from this term. That is why we keep the unity in the factor (1 + 2Nq) 
under the sum. As a result if = the third term in Eq. (|26|l is equal to the sum of the last two terms in Eq. 124(1 . 
The first two terms in Eq. 1(2611 result in the following renormalization of constants D and wq in Eq. (Q: 



D 



1 - ^W{T) 



UJQ ^ UJQ 



1 - ^V{T) 



(29) 



At the same time the last term in Eq. ((26(1 has another structure. It contains the isotropic correction to the spin- wave 
stiffness. As it is explained above, such a correction plays the predominant part when sin 0k = and k <C luq/D. The 
last term in Eq. ((26(1 contributes also to the spin-wave gap. 

Notably, thermal contributions to the isotropic correction to the spin-wave stiffness and to the gap depend on the 
demagnetizing factor Afz . This surprising feature is attributed to the long-range nature of the dipolar interaction. 

It should be noted that only the first two terms in Eq. 1(26(1 have been obtained in Ref.'^, where only the regime 
was studied in which the thermal fluctuations prevail. The origin of this discrepancy is in the following. First, it was 
discussed the special geometry of the sample in Ref.'^: a macroscopic slab of finite thickness d with square faces of 
length L — > oo; the magnetization is parallel to the surface. Demagnetizing factor is zero under these conditions. 
Second, to put it in our terminology, the anomalous Green's functions F and -F^ were discarded in all sums on the 
path leading to Eq. ((23(1 . As a result the terms in Eq. 1(23(1 containing Bq under sums were missed. In particular 
the last term under the first sum in Eq. ((23(1 was discarded that gives rise to the second term in the second brackets 
in Eq. I(26() . i.e., that leads to the spin-wave gap and to the isotropic correction to the spin-wave stiffness if Afz — 0. 
That is why only renormalization of constants D and loq given by Eqs. I(29() was obtained in Ref.'^ from 7i4-terms. 

We point out that summation over small q is important in Eqs. ((27(1 and 1(28(1 . Thus terms in Eq. ((26(1 containing 
W{T) and V(T) appears also in the case of BCCL and FCCL. The difference is in constants w and v only. 

It should be noted that there is also a renormalization of D and ujq similar to that given by Eq. I(29() in quantum 
regime. But corresponding terms in f2*^*^ {uj, k) are much smaller than those presented in Eq. II24I) and we omitted 
them. 



B. Loop diagram 

We turn now to the loop diagram shown in Fig.^c). It appears from Tia-terms ((10(1 in the Hamiltonian. As a result 
of simple but tedious calculations some details of which are presented in Appendix ^ we obtain for the contribution 
to Q(uj, k) from this diagram: 



17(3) (^,k) ^ -Dk^^ ain^ 29^U{Lu,k,T) - {De + ^ Sin 



5^ 



-(l + 2iVq), (30) 



C/(w,k,T) 



Re(_Bqg*_)^) + gq£;q_k - 6^ 

eq[(eq + eq-k)2 - (w + i<5)2] 



(1 + 2iVq) 



2[Re(BqB*_k) + EqEq-k + eqeq-k][eq - Cq-k] 



[(eq + eq-k)2 ~{lu + iSy 



Eq-k 



)2 - {io + tsy 



(iVq-k - iVq) 



(31) 



where we set fc = in the sum of the second term in Eq. ((3()() . Dimensionless function L/(w, k, T) contains the infrared 
singularity obtained previously by Rahman and Mills. At A;, a; — > small q give the main contribution to the sum in 
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Eq. H31() . As a result one leads to estimations 

Re[f/(w,k, T)] 

lm[U{uj,k, T)] 




(32) 
(33) 



C. Resultant expressions 



Using Eqs. (|24|l . H26|) and l|3U|) one can derive the resultant expression for ri(a;,k). Quantum and thermal regimes 
should be considered: T < (Duj^Y^^ and T > {Dlo^Y^^. 
T < (Dtjg)i/3. From Eqs. ^ and (jSOI) we obtain 



^2 c^3 



Y = 



5 J ^|Q^^-2^Q-: 



(Qqi' + (Qr)' 



>' q 



4eo 



0.012, 



(34) 
(35) 



where X is given by Eq. (|25|l . the first term in the brackets in Eq. (|34|l originates from fi'^'*) (w, k) and to the rest terms 
contribute both r^'^^^(a;,k) and ri^''^(ijj, k). The gap which square is given by the last term in Eq. 134|l is proportional 

to tjQ^^sin^k- The logarithmic divergence of [/(w,k, T) is screened by the gap and contribution from the first term 

in Eq. H30|l was omitted in Eq. H34|l being of the form sin^ 

It is easy to see that Q^^ and can be replaced by Q^^ in Eq. (|35|) . After this transformation Y appears to be 
equal to C given by Eq. (|21|l . This coincidence is not accidental. As is discussed in Introduction, it looks reasonable 
that dipolar anisotropy given by Eq. H20I) should be accompanied with a spin-wave gap. The relation between the 
dipolar anisotropy and terms in Eq. (|34(l containing Y can be shown in the following non-rigorous way. Let us discuss 
large spins and try to take into account the dipolar anisotropy H20(l phcnomenologically by adding to the microscopic 
Hamiltonian JSJl the following expression [cf. Eq. H20|) ]: 



C 



Wo 



4^3 J ^ 



[isn'isff + isrfisn' + isffisn' 



(36) 



This term after Dyson-Maleev transformation ijHJl gives the contribution CwQ/(2J)aj^ak to the bilinear part (|5J| of the 
Hamiltonian that in turn leads to the shift £'k -Ek + CijJ^/{2J). Using this renormalization of i?k, equality C = Y 
and Eq. for the spectrum one recovers corrections in Eq. ()34|l containing Y . This consideration does not work if 
5 < 3/2 because Eq. is a constant in this case. 
T > {Dulf/^. One has from Eqs. ^ and (HI: 



rj(w,k) = - {Dk^f '^W{T)~ Dk^uoUm^ ei,\W{T) + V{T)]~2NzV{T)~Y 



^sin^^k (2M.V{T) + Y^ 



(37) 



H28I) and (|35|l . respectively. It is seen from Eq. H37|l that the 
k if A4 7^ and Nz = 0, respectively. It screens the divergence of 



where W{T), V{T) and Y are given by Eqs. ^ 

gap is proportional to woT^/^sinf^k and cjq^^ sin " 
U{uj, k, T). When = the spin- wave gap has the same form as in the quantum regime. 

It should be noted that the second term in the second brackets of Eq. (|26(l and the sum in the second term of 
Eq. pnji are proportional to uj^ and TD{ujo/D)^/'^ at T <C V^qD and T > VtooD, respectively. At small T they 
result in terms proportional to the constant Y in Eq. H34|) . The temperature corrections from these terms which were 
large at T ^ ^/lJqD cancel each other and the terms in Eq. I|34() proportional to Y survive in the thermal regime. 



D. Face- and body-centered cubic lattices 



In the case of FCCL and BCCL arguments similar to those presented at the beginning of Sec. IIII Al leads us to 
the conclusion that the diagram shown in Fig. ^ a) is zero. The general expressions (gg and ^ for riW(w,k) 
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and n'-^^uj,k) are valid for FCCL and BCCL. The results for n{uj,k) differ slightly from Eqs. §^ and At 
T < (Dcjg)^/^ one obtains 

^2 C 3 

niuj,^) = -Dk^-j (^isin''6lk + X2sin3 6lkC0s6'kSin3</)k-F) + F-^ sin^ 6ik, (38) 
where we assume that z- and x- axes are directed along (1, —1, 1) and (1, 1, 0), respectively, for BCCL one has 

X2 = V2Xi = ^ — cos((z, +qy + q,) « 0.017, (39) 

q 



for FCCL we obtain 



^ cjoOl^:^ 24eq ~ • : \ ) 



q 



and constant Y is given by Eq. 13 5() where components of the dipolar tensor are taken now within the coordinate 
frame with z-axis directed along the body diagonal of the cube. Numerical calculation gives 

r 0.03 for BCCL, 

\ 0.003 for FCCL. ^ > 

We see comparing Eqs. (I34|) and (|38|l that the term proportional to the square invariant in the case of SCL is replaced 
by the term which can be expressed via triangular invariant Im([A:2: + iky]^)/k for BCCL and FCCL. This result is 
quite natural because the quantized axis in BCCL and FCCL, body diagonal of the cube, is triad axis. 

In thermal regime Eq. ^ is vahd for FCCL and BCCL, where W{T) and V{T) are given by Eqs. lf?7 | and if^ . 
respectively. The numerical constants w and v in these expressions are following: w — 0.09 and 0.03 for BCCL and 
FCCL, respectively; v = 0.23 and 0.12 for BCCL and FCCL, respectively. 

The relation between dipolar anisotropy (|20|) and terms containing Y in Eq. H38|l can be demonstrated (non- 
rigorously) as it is done above for SCL. After addition expression (|36ll to Hamiltonian (jSJ its bilinear part (O acquires 
new term —Cojq / {3J)al.a]f^ that in turn leads to the following renormalization: Ey^ i-^ — CuJq/{3J). Considering 
this renormalization, using Eq. H14|) for the spectrum and taking into account the equality 

C=-Iy (43) 
that follows from Eqs. H21|) and H35|l one recovers terms containing Y in Eq. H38|) . 



IV. DISCUSSIONS 



As is seen from Eqs. H34(l and H37|l . the spin- wave gap screens the logarithmic divergence of the first term in Eq. H30(l 
making it negligibly small. Let us discuss briefly further 1/5-contributions. Bearing in mind that imaginary part 
of U{Lu,'k,T) given by Eq. H31|) possesses power-law singularity i^JloqI DT ju) one could expect that divergences of 
further l/S'-contributions to the stiffness contain powers of ^Jun/DT/uj. The spin- wave gap screens such singularities 
and the corresponding terms turn into powers oiT /uoq. Thus, it is clear that at T ^ ojq the analysis of higher order 
l/S'-terms is not required. But such an analysis is necessary at T ^ loq that is out of the scope of the present paper. 

Let us turn to the problem of infrared divergence of longitudinal spin fluctuations obtained in Ref.^. The longitudinal 
spin Green's function is given by X||(a;,k) = i ^^'"'([^^(t), S'^(O)]). Within the first order of l/S* one has X|| (w, k) — 
AU{u), k, T)/SJ. It is clear from the above consideration that the spin-wave gap screens the divergence of U(uj, k, T). 
The analysis of further l/S'-terms is not required at T <C wo- However such an analysis is necessary at T 3> wo- 

Finally, we discuss the possibility to verify by inelastic neutron scattering theoretical results obtained above. Unfor- 
tunately to the best of our knowledge there is no suitable substance for the corresponding nowaday experiment. The 
most investigated compounds in which dipolar interaction plays important role (wq ~ D) are europium monochalco- 
genides EuO and EuS. But according to the analysis of Ref>i&, one-ion anisotropy is very large in these materials 
and seems to screen all the peculiarities discussed above. Seemingly, there is no such problem in chalcogenide spinel 
CdCr2Se4. According to estimations made in Refi^ the anisotropy field observed in this substance is of the order of 
dipolar anisotropy. At the same time loq/D ^ 10~^ in CdCr2Se4 and extremely small impulses transfer of the order 

of 10~^ A should be achieved in the corresponding neutron experiment. It is hardly possible nowadays. 



9 



V. CONCLUSION 

In the present paper an evaluation is carried out of the real part of first 1 / S'-corrections to the spin- wave spectrum 
of a cubic ferromagnet with dipolar forces. In accordance with previous results^i we obtain thermal renormalization 
(|29|l of constants D and ujq in the bare spectrum (Q. Besides, a number of previously unknown features are revealed. 
We observe terms which depend on azimuthal angle of the momentum k. They are proportional to the square 
invariant — fc^|/fc in the case of simple cubic lattice (SCL) and to the triangular invariant lT[i{[kx + ikyY)/k'^ in 
face-centered cubic lattice (FCCL) and body-centered cubic lattice (BCCL), where we assume that z-axis is directed 
along magnetization. This difference is accounted for the fact that dipolar interaction gives rise to the anisotropy: 
edges of a cube are easy directions in SCL and body diagonals of a cube are easy directions in BCCL and FCCLiSiS 
It is shown that there is a correction to the spin-wave spectrum proportional to k which does not depend on direction 
of k (isotropic correction). Thus the spectrum appears to be linear in k rather than quadratic [cf. Eq. (^] when 
k <^ loq/D and sin 0k = 0, where is the polar angle of k. In particular, we observe a spin- wave gap proportional to 
sin(?k- We demonstrate a relation between the dipolar anisotropy and the gap and isotropic term in the spectrum. 

Notably, thermal contribution from the Hartree-Fock diagram shown in Fig.^b) to the isotropic correction as well 
as to the spin-wave gap are proportional to the demagnetizing factor in the direction of the domain magnetization. 
Then the spin- wave excitations with wavelengths much smaller than the size of the domain appear to "feel" the 
domain shape. This nontrivial behavior is attributed to the long-range nature of the dipolar interaction. 

It is demonstrated that the gap screens singularities of perturbation corrections to the spin-wave stiffness^ and 
longitudinal dynamical spin susceptibilitji obtained before. We show that higher order l/S'-corrections to these 
quantities are small at T <C i^-'O; where uq is the characteristic dipolar energy given by Eq. However the analysis 
of entire perturbation series is still required to derive the spectrum and the susceptibility when T ^ wq- Such 
an analysis is out of the scope of the present paper. One could expect large renormalization of the longitudinal 
dynamical susceptibility at T ^ ujq as the nontrivial power-law dependence of X\\i^ 0) ~ (i/cj)'^'^^ was obtained 
experimentally at such temperatures 
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APPENDIX A: CALCULATION OF !^(^)(tj,k) 

We present in this appendix some details of calculation of f7'^^^(a;,k) that is a contribution to Q{uj,'k) given by 
Eq. H18|) from the loop diagram shown in Fig. ^c). This diagram appears from Tia-terms (|10|l in the Hamiltonian. 
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As a result of simple calculations we lead to quit a cumbersome expression: 

S 



Dk'^Lu^ sin^ 26'k 



qi:q2 



'«i«i-k'«-2ak-2) 



S 



sm t^k 

2 ^ 



S 



2m ^ 

qi,q2 



qi,q2 

^ ^ [(Miai_k-^^-2a2-k)^- (Afiak-i,M_24-2)- 
qi:q2 

sin26'ke""^''^ ^ [(ai4_k, M-2a2-k)c^ - (Miflk-i, a-24 
S 



k-2) 



■ a;wosin26'ke"^''^ ^ [(MiaLk' a-2aL2>" - (aiOi-k' ^-^-24-2)'^ 
qi,q2 



i:ifc2^osin26lke'"^''^ ^ [(Miai_k, a_2aL2>" + (ai«Lk. ^-^-24-2) 



4^ 



Dfc2cuosin20ke-*^''— ^ [(aioLk' A^-2«2-k> 



(Miflk-i, 0-2^-2) 



qi,q2 



Sk— XI (^^i«k-i,M_2a2-k> 



qi,q2 



k-2/'^' 



(Ala) 
(Alb) 
(Ale) 
(Aid) 
(Ale) 
(Alf) 
(Alg) 
(Alh) 



where 



denotes retarded Green's function, indexes q are dropped in expressions under sums and combination 

= aq(g^^ - iQin + alJQl^ + zQl^) (A2) 



is introduced. For the first glance expression IjAip looks awesome. Meanwhile it becomes simpler if one takes advantage 
of the properties of the combination Mq. 



(Mq, a- 



y 9<i Wo {Dq^ - uj) sin 26'q 



-Bq(Qq^ + ^Ql') + {Eq + ^)(Qq^ " ^Q^') 



2 {u + i5Y - e2 
<?<! Wo {Dq^ + w) sin20q 



(w + i(5)2 



(A3) 



(Mq,M_q)^ = 2i;q 



<3«i Wq Dq^ sin^ 2^q 
T (w + i(5)2 - el 



where we use expressions and {TJ for components of the dipolar tensor. One can demonstrate that terms (c)-(h) 
in Eq. IjAip are small compared to il(^'(w,k) given by Eqs. H24|l and H26|l . It is easy to see that summation over small 
momentum give the main contribution to these terms. Their evaluation using Eqs. (|A3() is a tedious but straightforward 

work. Terms (h) vanish as fc — > or 0k ^ 0,7r. As a result they are proportional to k'^uj^^^ sva^ 9]^ and oj^sm^ 9-\^/k 
at A: <C ^/bJo/D and k 3> yuJo/D, respectively. Sums in terms (f) and (g) vanish as fc ^ or 0k — > 0, 7r/2, tt. In the 
issue (f ) and (g) appear to be small compared to O'^^ (w, k) being proportional to k'^ujy^ sit? 20k and k^ujQ sir? 20k at 
k ^ \J Lo^jD and k 3> ^Jloq/ respectively. Sums in (d) and (e) are proportional to to and vanish as 0k — > 0, 7r/2, tt. 
Then, (d) and (e) can be also discarded being proportional to uj^ ku! a sir? 29\^ and uj'^loqSvcP' 29\^/k at fc ^ \Jijj{)/D 
and k ^ ^Jujo/D, respectively. Term (c) is negligible being proportional to w^Wq^^ and uJ^w^jk at /c ^ ^Juiq/D and 
k ^ y^uJo/D, respectively. The above evaluations are valid for T = 0. It can be demonstrated that temperature 
corrections to (c)-(h) are also negligible compared to il,^'^^uj,'k). 

As a result only terms (a) and (b) in Eq. (|Aip should be taken into account. They contribute to the first and 
the second terms in Eq. (ISOJ, respectively. This result is valid for all types of cubic lattices because summation over 
small momenta give the main contribution to (c)-(h). We point out that only normal self-energy parts (S and S) 
contribute to terms (a) and (b). 

Notice that term (a) is the only source of singularity in Eq. (|A1() . As is seen from Eqs. (IA3|) . the combination Mq 



11 



o 



a) 



b) 




FIG. 1: Diagrams of the first order of 1/ S. Diagrams a) and c) stem from three-magnons terms 1101 in the Hamiltonian whereas 
b) comes from four-magnons terms Ijll^ . 



is "soft": numerators of all the products containing Afq vanishes as g, w — > 0. 
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